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^•^Cu and ^®^Hg NMR shifts for an optimally doped and underdoped HgBa2Cu04+4 single crystal 
are reported, and the temperature dependence dictates a two-component description of the uniform 
spin susceptibility. The first component, associated with the pseudogap phenomenon in the NMR 
shifts, decreases already at room temperature and continues to drop as the temperature is lowered, 
without a drastic change at the transition temperature into the superconducting state. The second 
component is temperature independent above the superconducting transition temperature and van- 
ishes rapidly below it. \t is a substantial part of the total T dependent susceptibility measured at 
both nuclei. 



I. INTRODUCTION 

One of the most challenging questions in condensed 
matter theory is embodied in the physics of high- 
temperature superconducting cuprates and concerns 
their transition from a Mott antiferromagnetic insulat- 
ing state into a conductor and superconductor by remov- 
ing or adding electrons, i.e., by doping them. Clearly, 
understanding the electronic spin or magnetic spin sus- 
ceptibility xs of these materials is of great interest. In 
particular, it is important to know whether the uniform 
spin susceptibility can be understood in terms of a sin- 
gle electronic spin component that changes with dop- 
ing and temperature across the phase diagram. Nuclear 
magnetic resonance (NMR) measures xs locally through 
the hyperfine coupling that causes, e.g., a relative shift 
Ks,ri — Pi) • Xs of the nuclear resonance frequency, where 
p,, is the anisotropic hyperfine coupling constant^ and 
denotes the direction of the magnetic field with respect 
to the crystal axes. With NMR as a bulk probe, such 
measurements are reliable since they are insensitive to 
dilute impurities or magnetic inclusions. 

If all nuclei, n, of a given material experi- 
ence the same normalized, proportional shift change 
A"fCs,r)/"p,, = Axs in the temperature interval AT, this 
is evidence for a single-component response. On the other 
hand, if one finds experimentally that this shift change 
with temperature differs for different nuclei, this proves 
the failure of a single-component view. This important 
question was addressed with NMR measurements early 
on. And it was concluded, rather quickly, with exper- 
iments on YBa2Cu3 06.63^ and later YBa2Cu40g- that 
a single-component description of the cuprates is appro- 
priate. However, pioneering measurements of the total 
susceptibility indicated two-component behavior of dif- 
ferent cuprates-, and NMR relaxation studies did not 
agree with a single-component model^. 

Recently, some of us showed with extensive NMR 



shift measurements that Lai.85Sro.i5Cu04 does not 
show a single-component spin response^, and since 
high quality IIgBa2Cu04+5 single crystals had become 
available^!^, we address here the important question 
whether the two-component behavior is a mere peculiar- 
ity of Lai.85Sro.i5Cu04 or a more universal property of 
the cuprates. 

We have acquired a comprehensive set of NMR data 
for an underdoped and optimally doped sample^ , and we 
will show that both materials require a two-component 
description. Concurrent with the present work, some of 
us were involved with developing a new anvil cell NMR 
probe^*' for sensitive, high-pressure NMR experiments to 
enable, e.g., the investigation of the cuprates under pres- 
sures that had not be done before. Recent success was 
reported with ^^O NMR of YBa2Cu408 at pressures up 
to 63 kbar— . It was found, based in part on the results 
presented here, that the spin susceptibility in this system 
also requires a two-component description. 



II. EXPERIMENTAL 

The two crystals- used in this study were annealed to 
result in one optimally doped crystal with Tc = 97 K 
(mass 30.3 mg) that we label X97 and one underdoped 
crystal with Tc = 74 K (mass 3.3 mg) that we label X74^. 
^■^Cu and ^^^Hg NMR shifts as a function of temperature 
T were recorded in a magnetic field Bq = 11.75 T for 
two different orientations of the tetragonal single crystals 
with respect to the external magnetic field: with the field 
parallel (c [| Bq) and perpendicular (c_Li?o) to the crystal 
c-axis. NMR shifts were referenced using a fine metallic 
copper powder. For ^"^Cu we used the published value of 
the 63Cu Knight shift of 3820 ppm.^^ i99Hg shifts were 
referenced to (CH3)2Hg using the referencing procedure 
described by HarrisH. Below Tc, in particular for ^^^Hg, 
the sensitivity loss due to radio frequency penetration 
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FIG. 1. T dependencies of the total magnetic shifts ^^J^y ^(T) 
for both crystals, X74 (Tc = 74K ) and X97 (Tc = 97K); a) 
c II Bo, and b) c_LBo. The quadrupole interaction was deter- 
mined and its effect on the shift removed for c_LBo. 



FIG. 2. ''^-K'ii(r) versus ^^K^{T) with T as an imphcit 
parameter, for both single crystals, X74 and X97. The two 
arrows indicate Tc for both crystals. The two straight lines 
are linear fits to the data over the entire range for X97 and 
above Tc for X74. 



and slow nuclear relaxation limits data collection. Since 
the materials are in a mixed superconducting state below 
Tc, a macroscopic diamagnetic contribution Ky> that is 
not precisely known appears. Since K^) is isotope inde- 
pendent for a given sample orientation, it can be removed 
by taking shift differences, e.g., ^^Kr^ - ^^^i^T^. 



III. RESULTS AND DISCUSSION 

Fig. [T] shows the T dependence of the total magnetic 
Cu shifts ^^K,^ (the error bars are within the size of the 
symbols) that are the sum of orbital, if orb, and spin, Ks^ 
contributions. Note that this cuprate is known to show a 
T dependent Cu shift also for c || B[r^"— , in contrast to, 
e.g., YBa2Cu307_5 and La2-a;Srj;Cu04. The term due 
to the partly diamagnetic response in the mixed state, 
i^D, that we expect below Tc, can be neglected for the 
large Cu shifts, as we show below when we discuss the 
i99Hg NMR data. Note that the anisotropics of the shifts 
(compare the actual numbers in both panels in Fig.[T]) are 
due to the orbital (Van Vleck) susceptibility. We show 
the total magnetic shifts (verified with field dependent 
measurements) to keep the discussion transparent. 

In order to isolate the spin shift one usually assumes 
that -K^orb is largely T independent (as are the quadrupole 
interaction and the hyperfine coefficients), and that the 
T dependence of the shift is given by a T dependent spin 
shift that vanishes below Tc. This means xs(?^ — 0) 
due to spin singlet pairing, and thus the low-T shifts must 
be given by the orbital term. 

The fact that the shifts for X97 and X74 at the low- 
est T and given orientation are similar for each ori- 
entation in Fig. [T] shows that the orbital shift term 
does not depend significantly on doping (similar to, e.g.. 



La2-2:Sra:Cu04ii). At higher T the T dependence must 
then be due to the T dependent spin susceptibility. 

Let us assume momentarily that the shifts in Fig. [1] 
are caused by a single component spin susceptibility 
Xs{T). Then, all four curves in Fig. [T] are described 
by ^^Ks,^{xj,T) = 63p^^ . xs{xj,T) where ?? denotes 
the orientation of the crystal with respect to the mag- 
netic field and xi, X2 denote the doping. If we plot 
for a given sample the shifts measured for one orien- 
tation against that measured for the other orientation 
we expect straight lines with similar slopes. This can 
be seen by substituting xs{xj,T), e.g., ^^K^^\\{xj,T) — 
^\\/^^V±-^^Ks.x{xj,T). Such plots are shown in Fig. [2] 
(the actual data are the same as those presented in 
Fig. [T]). We observe indeed straight lines with similar 
slopes at higher T (larger shifts), however, there is an 
offset between them. Note that the similar slopes prove 
that the ratio of the hyperfine coefficients, and thus most 
likely the coefficients themselves, is independent on dop- 
ing (the same conclusion is reached when comparing Cu 
and Hg data, see below). For X97 the straight line con- 
tinues through Tc down to the lowest temperatures, while 
for X74 the measured data points move away from the 
straight line near Tc to meet those for the X97 material at 
the lowest T. That they meet at low T is expected from 
Fig. [U the orbital shifts are the same for both samples. 
Note that while we expect constant slopes, the constant 
high-T offset between the two lines, which disappears at 
low T, cannot be explained. 

Clearly, the observation in Fig. [5] that the changes of 
the X74 shifts with T are not proportional to each other 
at all T (here for different orientations of the field rather 
than different nuclei) demands a second spin susceptibil- 
ity with a different coupling to the nucleus. Thus, we are 
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led to write^ 

K^{xj,T) ■xi(xj,T) + qr, ■X2{xj,T), (1) 

where X2{T) is the second T dependent spin susceptibil- 
ity and the corresponding hyperfine coefficient. For a 
particular crystal {xj ) we then have for the data in Fig. [31 



KAT)='^-K^{T) 



qil 



• qi 



X2(T). (2) 



Then, X2(X97,T) must be rather small and X2(X74,T > 
Tc) = const., but vanishes below about Tc, and both sets 
of data in Fig. [2] meet at low T. 

The need for two shift components is already obvious 
in Fig. [T]where the shifts for both orientations of the X97 
crystal are drawn to scale in panels a) and b). Then, since 
the shifts in each panel are the same at low T the shifts 
for both orientations of the X74 crystal (lower curves in 
both panels) should look very similar on this scale, as 
well. This is clearly not the case and Fig. [1] tells us that 
the shifts cannot be understood with a single-component 
spin susceptibility. Instead, we need at least two com- 
ponents with different T dependences. We can also con- 
clude from the data in Fig. [T] that both components have 
to disappear at low T. 

We would like to point out that the two susceptibilities 
that couple with p^ ^ to the nucleus, cf. ([!]), may not 
be the suceptibilities of the actual spin components. For 
example, a second spin component that couples to the 
nucleus through the same p,, as the first one could be 
present in the shift data for the X97 crystal in Fig. [T] 
while not showing up in Fig. [2] However, we do infer 
from Fig. [5] that X2 is independent of T above T^. for 
the X74 crystal. We refer to a more detailed discussion 
further below. 

We now turn to ^^^Hg NMR. We show in Fig. [3] data 
for the X97 crystal: a pair of typical spectra for both 
orientations of the field (a), and the T dependencies of the 
corresponding shifts (b, c). We observe two well resolved 
i99jjg lijj^gg label them with a and /3. Note that the 
T dependence of the Hg shifts is not very different from 
that of Cu, but the shifts are smaller, as one expects 
from transferred hyperfine coupling to the more distant 
Hg atoms. Similar results were obtained for the X74 
crystal (data not shown) , which reveal that the a line is 
the same as for the X97 crystal at the lowest T where the 
spin shift is small, whereas the /3 line has a smaller shift 
and lower intensity. Therefore, the /? Hg site must be 
due to Hg atoms that are affected by doping (for further 
details we refer to^). Here we would like to focus on the 
X97 crystal in order to test the single vs. two component 
description for this system. 

Since we have two sharp, well resolved Hg resonances 
we can plot their shifts against each other, similar to what 
we did for ^^Cu, to inquire about the spin susceptibility. 
Note, however, that the orbital contributions to the shifts 
are very large compared to the T dependent shifts, see 
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FIG. 3. ^^^Hg NMR of the Tc = 97K single crystal (X97); 
a), typical spectra for both directions of the field show two 
distinct Hg sites; b) and c), T dependencies of the shifts for 
the lines shown above for c || Bo and c±Bo, respectively. 
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FIG. 4. X97 crystal: a) , b) T dependence of 
- 1.4 • i^^if,,!! and "^T^^j.x - 0.8 • "^if^.x, 
respectively. c) A^^'Ti'x = ^'^'^Ki3,i_ - ^^'^ K^^^ vs. 
A^''^Jf|l = ^'^^Kp t^\ — ^^^A'cii with r as an implicit parame- 
ter (see text). 
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both panels in Fig. |3l It is therefore most reUable to plot 
the shifts of both lines a and /3 for each orientation t] 
of the field against each other, i.e., ^^^X^^^ vs. ^'^'^Ka^-q- 
From such plots (see below) we find straight lines with 
fixed slopes for T > T^. Since K^^^ for T >Tcwe infer 
that there is only one T dependent spin susceptibility for 
T > Tc. We find for the slopes t\\ = 1.4 and r^ = 0.8 for 
c II Bq and C-LBq, respectively. In order to demonstrate 
this with higher graphical resolution we plot in Fig. |4] (a. 



b), ['''Kf,,,iT) 
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© to discuss Fig. g] (a, b). Since r^ 



Ka.jjiT)] vs. T. We can look at 
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horizontal lines for T > Tc in Fig.|4](a, b) are equivalent, 
cf . © , to constant slopes r^ above T^. in the correspond- 
ing shift-shift plots. Below Tc we observe deviations from 
the horizontal lines for both orientations, a clear indica- 
tion of another shift contribution. However, we cannot 
conclude from it with ([2|) on a second T dependent spin 
component since Kd could cause such a behavior, as well. 

In order to remove the uncertainty due to Kj) we now 
take the differences between the shifts of the two lines 



for each orientation. With A^^^K^ = ^^^if/3,r, - ^^'^i^a.r, 
we have with ([T]), here adopted for the ^^^Hg NMR spin 
shift. 
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XiiT) 



X2(r), (3) 



where A^'^-'p^ = P/3_,, — Pq.,,, etc. Analogous to 
([2) we can express one shift difference in terms of the 
other in order to remove one of the susceptibilities, and 
we have (dropping the mass number of the isotope for 
convenience) 



Apii 



Aqj 



Ap 



Ap 



X2. (4) 



This is shown in Fig. |4](c), however, there we plot the 
difference of the total experimental shifts between the /3 
and a line for c_LBo as a function of the same difference 
for c II Bq. Note that T independent orbital terms cause 
a constant offset on both axes. For a single T dependent 
spin component, i.e., X2 = in a straight line with 
a fixed slope is expected. This is clearly not the case. 
We observe that the shift difference for c || Bq drops sig- 
nificantly already at high T as the T is lowered, while 
that for c_LBo is constant at high T, but drops only be- 
low about Tc, a behavior that cannot be explained by 
a single T dependent spin susceptibility. From (|4]) we 
conclude that we need a second susceptibility that is T 
independent above Tc, but changes below Tc. This inde- 
pendent experimental observation proves that a second 
component is necessary also for the X97 crystal even if 
it appears to be small. While this is in agreement with 
what we derived from comparing the Cu shifts for X97 
and X74, it was not obvious from the Cu shifts for the 
X97 crystal alone. 

Somewhat unusual plots are presented in Fig. [Sj In 
the main panel we see the Cu shift with c || Bp for one 
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In order to remove a possible diamagnetic term Ko 

the Hg shifts have been subtracted from the Cu 

shifts in the inset: [''^-K'n X97(T) - "^i^ii X97(r)] vs. 

[«3if||,X74(r)-"'if||..,X74(r)i. 



crystal (X97, vertical axis) plotted against the Cu shift 
with c II Bq of the other crystal (X74, horizontal axis). 
Surprisingly, one finds a straight line above the X97 Tc. 
This says that, coming from high T (large shifts), the 
shifts in Fig. [T] for both crystals fall in such a way that 
leaves the ratio of the slopes constant. At the Tc of the 
X97 crystal its shift drops rapidly and below 70 K we 
find a behavior similar to the one observed at higher 
T. The fact that they arrive at similar shift values for 
the lowest T has been noted in Fig. [T] and is expected 
for a doping independent orbital shift. Since we know 
that X2(r > Tc) — const, we conclude from Fig. [5] for 
T > Tc that since [Ai4rx97] / [Ai^rx74] ~ 1.4 we must have 
Axi(X97)/Axi(X74) w 1.4 for aU T intervals. Thus, 
Xi(X97) = 1.4 • xi(X74) + c, where c is a constant that 
must have vanished below 70K in Fig. [5l This says that 
Xi consists of a component that is T dependent at all T 
and changes proportionally with doping, plus a second 
component with a behavior similar to X2- 

Before we continue with the discussion we would like 
to prove that the diamagnetic shift Kd is indeed not im- 
portant when comparing Cu shifts. We show in the inset 
of Fig. [S] a plot similar to that in the main panel, but with 
the Hg shifts subtracted from the corresponding Cu shifts 
so that K^) is eliminated on both axes. As expected, we 
observe again a straight line above Tc and since the Hg 
shifts are small it has a similar slope. The drop just be- 
low Tc remains and this rules out the importance of Kj:, 
on the Cu shift scale. 

We plot in Fig. [S] the Hg shifts (total experimental 
shifts) of the a lines against those of Cu for both crystals 
and c II Bq. As expected, the slopes at higher T are 
approximately straight lines and similar (note that larger 
errors are expected since the Cu shifts are much bigger 
than those of Hg). We find that ^^^P||,a/^^P|| ~ 0.12. 
For the X97 crystal we observe a steeper drop near Tc, 
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FIG. 6. X74 and X97 crystals: ^^^K« 
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crystals with T as an implicit parameter (see text). 



which levels off to the slope above Tc before it finally 
drops much more rapidly with lowering T. Since the 
first drop below Tc for the X97 crystal is accompanied 
by a large shift change for Cu, it is not dominated by 
Ky). The second and steep drop for both crystals at low 
T must be mostly due to Ky^ since its slope is close to 
1 {Ky) is independent of the nucleus) . With this we can 
estimate the maximum diamagnetic shift for c || Bq to 
be Ai^D < 0.015%. Thus, the contributions to the T 
dependent shift below Tc from Kd and are of similar 
size for ^^^Hg (but not for ^"^Cu). 

Since we need two spin component (we call them A 
and B) for the explanation of our data, we would like 
to point out that this demands in general^* three spin 
susceptibilities^: xaa,Xbb, and xab since there can be 
a coupling of the two spin components A and B. Accord- 
ingly, we have xi = Xaa + Xab, and X2 = Xab + Xbb- 
The total spin susceptibility is then given by xs = 
XAA + 2xab + Xbb- For the NMR spin shifts it follows, 

Kn,x = Pr, (XAA,x + XAB,x)+q,, {xbb.x + XAB,x) • (5) 

We know that (xbb,x + Xab,x) must be T independent 
at higher T. We now assume that both component xbb.x 
and xab,x indeed contribute to the shifts. Since it is 
highly unlikely that both components are T dependent, 
but their sum is not, we conclude that both must be T 
independent above Tc. Then, xaa must be the suscepti- 
bility that is T dependent already above Tc. 



We deduce from Fig. [2] that P||/p± ~ 0.40 and with 



Cx = [q|| - (P||/p_l) • q±] [xab,x + Xbb,x] 



(6) 



that Cx97 ~ 0%, Cx74(T > Tc) w -0.040%, 
and Cx74(T < Tc) « 0%. Thus, at high T 

[xab,X74 + Xbb,X74] ~ const., but it vanishes rapidly 
below Tc. We also find xab,X97 ~ -Xbb,X97, and if 
one of the components is sizable this relation holds 



approximately also below Tc. We now understand why 
the influence of the second spin component is rather 
weak in the plots in Figs. [2]and|4] (c) for X97. 

Note, that we discovered with Fig. [S] that 
AKxgr / AKx74 ~ 1.4 at higher T. A simple explanation 
follows by assuming that xaa,X97 ~ 1.4 • xaa,X74 and 
that the integration constant is just due to xab,x- We 
then have with ^ for the curve in the main panel. 
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^^r;,X97 = Kjj^XTA 

+ ^^Pr) ■ XAB,X97 



(7) 



+ ^^q>) • [XAB,X97 + XBB,X97] 
~ S {i^^Pv + ^^qi)) XAB.X74 + '^^qn ■ XBB.X74} ■ 

We note that the terms in the last row stem from the 
X74 cyrstal and will not contribute to the drop in Fig. [5] 
near Tc « 97K. In addition, since the X74 shift for c || Bq 
in Fig. [1] does not change significantly below Tc we know 
that the term in the last row can be ignored at lower T, 
as well. Similarly, the term in the third row was found 
to be rather small for the X97 crystal (x2,X97 ~ 0). We 
are left with the term in the second row, and since xab 
vanishes rapidly at Tc « 97K the observed drop in Fig. [5] 
is expected. We conclude that this drop is mainly given 
by 63p|| . AXAB,X97 ~ +0.05%. 

We can also inspect Fig. [5] again. For the X97 crystal 
the change of the slope at Tc marks the change in xab 
and Xbb until A'd dominates at lower T. On the horizon- 
tal axes this is given by A^^K^xgr = *^^P|| ' Axab,X97 ~ 
H-0.05%, as concluded from Fig. [S] The corresponding 
change in shift for Hg of about +0.01% is somewhat big- 
ger due to Ku C^^P\\/^^P\\ ~ 0.12). For the X74 crystal 
we do not expect a drop near Tc for Cu (cf. discussion 
above), and for Hg it is mostly given by Kb- (It appears 
that the onset of the change in slope is slightly above Tc 
outside the estimated error in the T measurement). 

On general grounds we expect xaa to be positive and 
therefore it follows that ^"^p^ > for both orientations. 
Then, we have for T > Tc that xab,X97 > 0, and since 
XAB,X97 ~ -XBB,X97, it foUow that xbb,X97 < above 
Tc for this material. For the X74 cyrstal we notice, e.g., 
from Fig. [U that the constant shift term above Tc must 
be small for c || Bq, but larger for c+Bq. 

To conclude, based on a thorough NMR study of the 
two HgBa2Cu04+5 crystals^ii^ we identified a 2nd shift 
component that is independent of T above Tc and van- 
ishes rapidly below Tc. Since it influences the Cu and 
Hg nuclei by similar relative values, the 2nd component 
must be due to a spin susceptibility, in accord with what 
some of us proposed earlier for Lai.85Sro.i5Cu04^, and 
very recently for YBa2Cu40g^^. Our results can be ex- 
plained consistently with two spin components that are 
coupled to each other and thus have the three suscepti- 
bilities XAA, XAB, Xbb- In such a scenario, cf. ([5]), we 
find that only one of them, xaa is T dependent above 
Tc, while the other two vanish below Tc. Furthermore, we 
find that XAA.X97 = l-4-XAA,X74, a relation that has been 
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suggested by some of us earlier^^ and that some of us ap- 
plied successfully to understand the changes in the ^^O 
NMR spin shift of YBa2Cu408 with pressureil. Thus, 
Xaa{T) is responsible for the pseudogap feature observed 
in NMR and it changes proportionally (increases) with 
increasing doping for an underdoped cuprate, but de- 
creases with increasing pressure^^. (We believe that 
XAA(a;2)/xAA(a;i) = X2/X1 where Xj is the local hole 
doping per Cu02 in the plane, so that xaa = p{T) ■ x, 
seei^). The part of the shift that is constant above 
was found to increase dramatically with pressureii and 
we believe that it mus also increase with doping. How- 
ever, more experiments are necessary to explore the two 
components across the phase diagram and as a function 
of pressure. Nevertheless, there can be no doubt that a 
two-component scenario is generic to the cuprates. This 
also means that the hitherto adopted single-component 
interpretation of NMR data can no longer be upheld. 

We can only guess the physical origin of the two 
spin components. Barzykin and Pines in an extensive 



review-- have put forward a model. Although the T de- 
pendence of the first component reminds one of coupled 
Heisenberg spins^, we note that the pseudogap regime of 
the cuprates has recently been shown to be associated 
with unconventional magnetism^li^^. Whereas the 2nd 
component exhibits the behavior of a Fermi liquid that 
becomes superconducting at Tc, a T constant behavior 
at higher T is also expected for a Fermi glass^. 
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